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Abstract 


Two  methods  are  developed  to  average  V  and  /  Stokes  parameters  over 
a  circular  sensor  aperture  collecting  light  scattered  from  an  optically  active 
sphere.  One  method  uses  a  two-dimensional  numerical  integration  of  Bohren’s 
theory  in  scalar  form,  which  explicitly  shows  the  connection  with  Mueller  ma¬ 
trix  elements.  The  second  method  uses  expansions  of  vector  spherical  har¬ 
monics  in  Bohren’s  theory  to  integrate  over  apertures  of  any  size  without 
convergence  checks  since  numerical  integration  is  avoided.  Equations  simpli¬ 
fying  the  analytical  results  for  dvr  scattering  are  obtained.  Sample  computa¬ 
tions  of  average  Stokes  parameters  are  performed  with  both  methods,  and 
agreement  to  six  decimal  places  is  obtained. 
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1.  Introduction 


The  presence  of  chiral  molecules  in  biological  materials  often  makes  these 
materials  optically  active:  i.e.,  characterized  by  different  complex  indices  of 
refraction  for  right  and  left  circularly  polarized  light.  A  material  with  dif¬ 
ferences  in  the  real  part  of  the  refractive  indices  is  referred  to  as  circularly 
birefringent,  and  a  material  with  differences  in  the  imaginary  part  is  referred 
to  as  circularly  dichroic  [1].  Many  optically  active  materials  exhibit  both 
circular  birefringence  (i.e.,  optical  rotation)  and  circular  dichroism. 

Early  laser  light  scattering  research  in  the  area  of  optically  active  particles 
includes  the  experimental  measurements  of  Holzwarth  et  al  [2]  in  the  1970s, 
who  used  an  application  of  Mie  theory  developed  by  Gordon  [3]  to  simu¬ 
late  scattering  from  spherical  microorganisms.  In  the  1990s,  measurements 
of  Mueller  matrix  elements  for  light  scattering  from  nonspherical  microor¬ 
ganisms  have  been  reported  by  Lofftus  et  al  [4]  and  by  Bronk  et  al  [5],  who 
also  performed  realistic  light  scattering  calculations.  These  results  provide 
evidence  that  laser-based  biosensors  can  detect  and  characterize  biological 
organisms  in  water  or  the  atmosphere  if  experimental  measurements  are  sup¬ 
ported  by  realistic  simulations. 

Exact  theories  of  light  scattering  from  optically  active  particles  were  hrst 
given  by  Bohren,  who  gave  correct  solutions  for  spheres  [1,6],  spherical  shells 
[7],  and  cylinders  [8].  The  theory  for  layered  cylinders  was  later  given  by 
Kluskens  and  Newman  [9],  and  the  theory  for  spheroids  was  given  most  re¬ 
cently  by  Cooray  and  Ciric  [10].  These  theories  provide  algorithms  for  scat¬ 
tering  in  a  single  direction.  Two  Stokes  parameters,  referred  to  here  as 
and  can  be  calculated  as  functions  of  the  scattering  direction,  and  the 

degree  of  circular  polarization  is  the  ratio  of  these  Stokes  parameters. 

In  earlier  work  [11],  two  of  us  (Rosen  and  Pendleton)  applied  Bohren’s  scat¬ 
tering  theory  in  the  hrst  simulation  of  scattering  from  optically  active  spheres 
into  sensor  apertures.  The  experimental  conhgurations  considered  in  that 
work  [11]  are  characterized  by  the  presence  of  a  sensor  with  a  circular  aper¬ 
ture.  The  assumption  is  that  all  light  scattered  into  the  solid  angle  inter¬ 
cepted  by  the  circular  aperture  is  focused  onto  the  sensor.  This  means  that 
we  obtain  realistic  simulations  of  the  measured  signals  not  by  computing  the 
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Stokes  parameters  for  a  single  direction  but  by  obtaining  the  average  of  each 
Stokes  parameter  over  the  aperture  solid  angle.  Our  earlier  paper  [11]  very 
briefly  outlined  the  theory  but  emphasized  a  discussion  of  computational  re¬ 
sults,  omitting  any  discussion  of  the  algorithm  used  for  the  computations. 
To  provide  a  better  basis  for  future  computations  of  this  type,  we  discuss 
the  theory  of  two  relevant  algorithms  here.  It  is  assumed  that  most  readers 
of  this  article  will  be  interested  in  developing  computer  programs  for  simu¬ 
lation  of  experiments,  so  we  have  taken  care  to  give  results  allowing  efficient 
computation. 

The  earlier  computations  [11]  were  performed  with  a  numerical  integration 
over  the  aperture  solid  angle.  We  developed  this  method  by  adapting  a  two- 
dimensional  Legendre-Gaussian  numerical  integration  algorithm  originally 
designed  by  Pierce  [12]  for  integration  over  a  circle.  The  advantage  of  this 
approach  over  less  systematic  techniques  is  that  it  gives  very  precise  and 
efficient  results,  because  the  location  of  the  grid  points  in  the  aperture  solid 
angle  is  optimized.  A  succinct  description  of  the  method  is  given  in  section 
3.  A  computer  program  based  on  this  method  was  developed  and  run  on 
a  personal  computer,  giving  the  computational  results  reported  earlier  [11]. 
This  numerical  integration  method  is  computationally  simple,  but  the  results 
must  be  checked  for  convergence;  i.e.,  the  computation  must  be  repeated  with 
an  increasing  number  of  points  in  the  computational  grid  until  convergence  is 
obtained.  Also,  this  method  cannot  treat  aperture  solid  angles  approaching 
27r  steradians. 

To  overcome  these  limitations  of  the  numerical  integration  method,  we  de¬ 
veloped  a  second  method  to  treat  any  aperture  solid  angle  without  conver¬ 
gence  checks.  This  method,  which  does  not  use  any  numerical  integration,  is 
referred  to  here  as  the  analytical  method.  Explicit  expressions  for  the  two  av¬ 
erage  Stokes  parameters  are  given  that  describe  the  scattering  of  a  wave  from 
an  optically  active  sphere  into  an  aperture  located  in  any  specihed  direction 
and  intercepting  a  solid  angle  of  any  size. 

The  analytical  method  is  applied  here  to  optically  active  spheres  for  the  first 
time.  The  averaging  of  is  an  easy  extension  of  an  approach  developed 

for  Mie  scattering  into  apertures  by  Chylek  [13],  Wiscombe  and  Chylek  [14], 
Chu  and  Robinson  [15],  Pendleton  [16,17],  and  Son,  Farmer,  and  Giel  [18,19]. 
The  averaging  of  with  this  method  is  new  and  is  the  aspect  of  this 

report  most  likely  to  appeal  to  theorists. 

The  theory  given  in  these  last  five  references  used  the  notation  of  Jackson 
[20].  We  originally  developed  the  theory  given  here  using  Jackson’s  notation 
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but  rewrote  it  using  notation  that  is  more  standard  in  light  scattering  [1] 
to  interface  with  Bohren’s  theory  for  spheres.  This  choice  of  notation  should 
allow  easy  application  to  other  particle  shapes  for  which  scattering  solutions 
are  known. 

Both  the  numerical  integration  and  analytical  methods  can  treat  an  incident 
wave  that  is  linearly  polarized,  randomly  polarized,  or  unpolarized.  Mathe¬ 
matically,  a  randomly  polarized  wave  is  equivalent  to  an  unpolarized  wave, 
so  these  results  are  referred  to  here  as  unpolarized. 

The  objectives  of  the  numerical  integration  and  analytical  methods  are  the 
same,  but  the  approaches  are  so  different  that  it  is  surprising  that  the  two 
methods  give  results  that  agree  to  better  than  six  decimal  places.  This  agree¬ 
ment  allows  one  to  have  conhdence  in  the  precision  of  these  methods,  and 
the  availability  of  two  methods  allows  one  to  apply  the  more  appropriate  to 
a  particular  situation. 

A  computer  program  based  on  the  analytical  method  is  more  computation¬ 
ally  intensive.  Such  a  computer  program  is  more  difficult  to  write  and  requires 
more  computer  time  for  small  apertures  than  would  a  program  for  the  nu¬ 
merical  integration  method.  The  advantage  of  the  analytical  method  is  that 
it  can  treat  any  aperture  size  without  convergence  checks. 

It  is  hoped  that  application  of  these  methods  will  stimulate  the  develop¬ 
ment  of  realistic  simulations  and  new  experimental  techniques  based  on  light 
scattering. 
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2.  Basic  Definitions 


2.1  Nondimensional  Stokes  Parameters  and  the  Effective  Degree 
of  Circular  Polarization 


For  computational  convenience,  we  define  nondimensional  Stokes  parameters 

'^Stokes  jStokes 


-^Stokes  —  ^2^2  {Ir  -  h)  ^ 
line 

j Stokes  _  ^^2^2  i^R  S~  ^l) 

line 


(la) 

(lb) 


where  Ir  and  II  are  the  irradiances  of  the  right  and  left  circularly  polar¬ 
ized  components  of  the  scattered  light,  line  is  the  irradiance  of  the  incident 
electromagnetic  wave,  and  k  =  2'it/\  is  the  wave  number  in  the  medium  sur¬ 
rounding  the  sphere  where  the  wavelength  is  A.  In  the  far  held  where  Ir  and 
II  are  inversely  proportional  to  these  Stokes  parameters  are  functions 
of  scattering  direction  only  and  are  not  functions  of  r,  the  distance  from  the 
sphere’s  center  to  the  point  where  the  Stokes  parameters  are  measured. 

We  obtain  the  average  Stokes  parameters,  i^jstokes^^  averag¬ 

ing  the  Stokes  parameters  over  the  sensor  aperture  solid  angle  Ahl  with  the 
expression 

and  the  effective  degree  of  circular  polarization  for  scattering  into  the  solid 
angle  Ahl  is  then  dehned  to  be  the  ratio  of  these  average  Stokes  parameters. 
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This  method  requires  that  two  electronic  signals  be  measured.  One  of  these 
signals  is  proportional  to  the  radiant  power  =  {IincAQ/k^){V^^°^^^)  and 


the  other  is  proportional  to  j ■  If  the  electronic  re¬ 

sponse  function  of  the  system  is  dehned  to  be  the  ratio  of  these  two  electronic 
signals,  then  the  electronic  response  function  is  independent  of  line  and  turns 
out  to  be  equal  to  the  effective  degree  of  circular  polarization. 

In  bulk  material  irradiated  by  a  linearly  polarized  beam,  Ir  and  Jl  of  the 
transmitted  wave  are  different  only  if  the  material  is  circularly  dichroic.  The 
idea  motivating  this  research  was  to  extend  this  concept  from  bulk  to  a  single 
particle  scattering  into  a  sensor  aperture  so  that  a  nonzero  value  of 
would  indicate  the  presence  of  an  optically  active  particle.  In  our  earlier  work 
[11],  the  occurrence  of  a  nonzero  value  of  p'^ffecuve  optically  inactive 

particle  was  dehned  to  be  a  “false”  signal  of  optical  activity.  We  claimed  [11] 
that  no  false  signals  could  occur  for  some  aperture  locations  and  identihed 
these  aperture  locations.  Here,  we  prove  this  claim  for  both  the  numerical 
and  analytical  methods. 

2.2  The  Coordinate  Systems  and  Euler  Angles 

The  development  of  theory  for  scattering  into  arbitrarily  located  apertures 
is  facilitated  by  the  introduction  of  an  aperture  coordinate  system  located 
with  respect  to  the  beam  system  by  Euler  angles. 

The  convention  in  Mie  theory,  followed  by  Bohren  in  his  theory  of  optical 
activity,  is  that  an  [x,  y,  z)  Cartesian  coordinate  system  is  located  with  the 
origin  at  the  center  of  the  sphere,  and  the  ;2-axis  extends  in  the  direction 
of  incident  wave  propagation.  If  the  incident  wave  is  linearly  polarized,  the 
x-axis  extends  in  the  direction  of  the  incident  electrical  held.  The  spherical 
system  associated  with  the  {x,  y,  z)  Cartesian  system  is  the  (r,  0,  0)  system 
with  basis  unit  vectors  (f,  0,  0).  To  perform  the  integration  over  the  aperture 
solid  angle  AH,  we  introduce  a  second  coordinate  system  {x^,y^,z"^)  with 
the  same  origin,  but  with  orientation  determined  by  the  center  of  the  aperture 
solid  angle  such  that  the  z^-axis  extends  through  the  center  of  the  aperture. 
In  the  angular  momentum  literature,  the  convention  is  to  use  superscript 
primes  on  the  second  coordinate  system.  Here,  the  primes  are  replaced  by 
A  (for  aperture).  The  two  systems  have  the  same  origin,  and  the  relative 
orientations  are  dehned  by  Euler  angles  a,  (3,  and  7  as  commonly  used  in  the 
quantum  theory  of  angular  momentum  [21,22].  The  connection  between  the 
beam  coordinates  {x,y,z)  and  the  aperture  coordinates  {x^,y^,z^)  is  given 
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by  the  matrix  equation 


where 


( 

yA 

=  R{a(3'y) 

y 

\z^} 

U/ 

Rz{a) 


RM 


cos  a  sin  a  0 
—  sin  a  cos  a  0  ; 


V  0 

0 

ly 

^  cos  0 

0  -sin0\ 

0 

1 

0 

\  sin  0 

0  cos  0  / 

(3a) 


(3b) 


(3c) 


(3d) 


The  matrix  Rzi'j)  is  obtained  by  replacement  of  a  by  7  in  equation  (3c).  The 
Euler  angles  a  and  (3  are  the  0  and  9  angles  describing  the  direction  of  the  2:"^- 
axis  (through  the  center  of  the  aperture)  as  measured  in  the  (x,  y,  z)  system. 
For  a  circular  aperture,  the  light  scattered  into  the  aperture  is  independent  of 
the  third  Euler  angle  7,  which  is  a  rotation  about  the  A  superscript 

A  means  that  the  variables  (r,  9, 0)  and  basis  vectors  (f,  9,  0)  of  the  spherical 
beam  system  are  to  be  replaced  by  the  variables  (r,  9^,  cj)^)  and  basis  vectors 

{r,9  ,0  )  of  the  aperture  system.  For  a  circular  aperture  centered  on  the  2:"^- 
axis,  setting  dil  =  sin  6*^  d9^  d(f)^  with  0  <  0^  <  27r  and  6*^^  <  9^  <  9^^^ 
gives  the  integration  over  the  aperture,  where  9^^^  dehnes  the  outside  edge  of 
the  aperture  and  9^^^  defines  the  inside  edge  of  the  aperture.  The  solid  angle 
with  9^  <  is  covered  by  a  concentric  “beam  stop”  or  “beam  block” 
so  that  An  =  27r(cos6'An  —  cos9^^^).  As  the  name  suggests,  a  beam  stop 
is  usually  used  to  stop  unscattered  light  of  the  incident  laser  beam  from 
irradiating  the  sensor  when  the  sensor  aperture  extends  into  the  beam. 
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3.  Numerical  Integration  of  Stokes  Parameters 


In  the  “far-field  approximation,”  the  scattered  held  is  approximately  trans¬ 
verse  and  may  be  written  as  E  =  EqO  +  E^cf).  The  irradiances  of  the  right 
and  left  circularly  polarized  components  of  the  scattered  light  are  obtained 
with  the  relations  In/hnc  =  Er  ■  E*J{EincY  and  hllinc  =  Er  ■  El/{EincY, 
where  Er  =  {Eg  +  iE^)/  Er  =  {Eg  -  iE^)/^/2,  and  Einc  is  the  magnitude 
of  the  incident  electrical  held.  The  Stokes  parameters  dehned  in  equation  (1) 
can  then  be  expressed  as 


yStokes  p  2^^^‘^^sE^) 

{E^ncY  ’ 

(4a) 

jStokes  7  2  2  [  -^0  +  -^</>  ] 

{E^ncY  ■ 

(4b) 

3.1  Linearly  Polarized  Incident  Wave 


Using  the  exp(— icnt)  convention  for  the  assumed  monochromatic  time  de¬ 
pendence,  we  can  represent  the  linearly  polarized  incident  wave  described 
above  as  =  xEincexp{ikz),  and  Bohren  [1]  obtained  the  solution  for  the 
electrical  held  scattered  from  an  optically  active  sphere  as  a  function  of  the 
vector  spherical  harmonics  and  The  superscript  (3)  on  these 

vector  spherical  harmonics  indicates  that  the  radial  function  is  hY\kr),  a 
spherical  Hankel  function. 


Bohren  expressed  his  solution  for  the  scattered  electrical  held  as 


E  =  y 


,.„(2n  +  l) 


^  n{n  +  1)  L 


[m„)N 


(3)  , 

eln  ' 


-idnj'N 


_l_  p  -I- 

oln  W  Cniyi-eln  W 


-bn)M, 


(3) 

oln 


(5) 


The  coefficients  and  bn  reduce  to  the  Mie  coefficients  with  the  same  names 
when  there  is  no  optical  activity.  For  an  optically  active  sphere,  dn  =  —Cn, 
and  when  there  is  no  optical  activity,  dn  =  — =  0. 
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In  the  far-field  approximation  where  h\^\kr)  — >■  the  compo¬ 

nents  of  the  scattered  held  E  may  be  expressed  in  standard  form  as 


/  Ee\_  exp  ikr  (  S2  83^  (  cos  0  \ 
[-E<j))~  -^kr  [s,  sj  Vsin0; 


(6a) 


where  the  2x2  matrix  is  the  amplitude  scattering  matrix  with  elements 
Sk  {k  =  1  to  4)  [1,23].  For  an  optically  active  sphere,  Bohren  identihed  the 
scattering  matrix  elements  as 


^1  = 
^2  = 
^3  = 
,S4  = 


E 

n=l 

00 

E 

n=l 

00 

E 

n=l 

00 

E 

n=l 


2n  +  1 
n{n  -|-  1) 

2n  -I-  1 
n{n  +  1) 

2n  -I-  1 
n{n  +  1) 

2n  -I-  1 
n{n  +  1) 


4“  bfiTn)^ 

{^OjnTn 

{CnT^n 


(6b) 

(6c) 

(6d) 

(6e) 


where  7r„  and  are  functions  of  cos  9  that  are  dehned  in  Mie  theory  [1,23] 
(see  app  A).  Since  dn  =  — for  optically  active  spheres,  it  is  found  that 
5*4  =  —S3.  For  spheres  that  are  not  optically  active.  Si  and  S2  become  the 


functions  of  the  same  name  in  Mie  theory  and  5*4  =  —  S'3  =  0.  Substituting 
equation  (6a)  into  equation  (4a)  gives 

y Stokes  ^  cos(20)  -I-  543  sin(20),  (7a) 

where  S'41,  5*42,  and  5*43  are  Mueller  matrices  dehned  as 

,S4i=Im(,Si,S3*-,S2^:),  (7b) 

,S42  =  Im(-,Si,S3*-,S2^:),  (7c) 

,S43  =  Im(,Si,S;-,S3A:).  (7d) 


The  approach  used  earlier  [11]  was  to  see  if  measurements  of  ^  0 

could  be  used  to  indicate  the  presence  of  an  optically  active  particle.  Then,  for 
spheres  with  no  optical  activity,  it  would  be  required  that  _  q_ 


optically  inactive  spheres,  S'41  =  5*42  =  0  but  S'43  7^  0  in  general.  However, 
we  show  below  that  if  a  =  0  or  7r/2,  then  (S'43  sin(2(;/)))  =  0.  The  result 
is  that  =  0  for  optically  inactive  spheres  when  the  aperture  is 

centered  in  the  x-z  plane  (a  =  0)  or  the  y-z  plane  (a  =  7r/2).  These  are  the 
preferred  locations  to  center  the  aperture  if  it  is  considered  desirable  to  have 
(yStokes\^  _  g  ^  Sphere  with  no  optical  activity,  i.e.,  no  false  signal. 


Substituting  equation  (6a)  into  equation  (4b)  gives 

jStokes  ^  ^  ^Qg^20)  +  ,Si3  Sin(20),  (8a) 

where  S*!!,  S'12,  and  S'13  are  Mueller  matrix  elements  dehned  as 

S*!!  =  -(|S'ip  +  |S'2p  +  |S'3p  +  |S'4p),  (8b) 

5i2  =  + i-^2r  -  + i-^4r),  (sc) 

^13  =  Re(^l^:+52^3*).  (8d) 


It  is  necessary  to  compute  cos{2(j))  and  sin(20)  for  each  point  on  the  grid. 
We  do  this  by  rewriting  equation  (3a)  as 


(vA 

^  sin  6  cos  (j)  ^ 

^  sin  9^  cos  (j)^  ^ 

V2 

= 

sin  6  sin  0 

=  [R{af3-f)]  ^ 

sin  9^  sin  (f)^ 

cos  9  y 

cos  9^  ) 

and  recognizing  that  cos(20)  =  [(ui)^  —  ('y2)^]/[('i'i)^  +  ('^'2)^]-  The  other  re¬ 
quired  expression  is  sin(2(;/))  =  2viV2l[{viY  -|-  (^2)^]-  This  expression  can  be 
used  with  equation  (8)  to  show  analytically  that  if  a  =  0  or  7r/2,  then 
ly Stokes^  _  (5'^3  sin(2(;/)))  =  0  for  a  sphere  with  no  optical  activity. 

Equation  (2a)  will  now  be  adapted  for  numerical  integration.  A  differential 
element  of  area  on  the  flat  circular  aperture  surface  is  dA^  =  dp^  d(f)^  = 
{z^YdVt/  cos^  9^  so  that 

where  F  =  f{cos9^)^.  The  Legendre-Gaussian  numerical  integration  scheme 
of  Pierce  gives  the  algorithm  for  integrating  some  function  F  over  a  planar 
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annulus  (i.e.,  the  area  between  two  concentric  circles),  and  the  solution  to 
these  integrals  may  be  expressed  as 

1 1  dA^F  =  E  E  (11) 

*  j 

where  the  area  of  the  annulus  is  =  7r[pEx  ~  Pmin]-  The  Wij  are  weighting 
numbers  normalized  so  that  J2iJ2jWij  =  1,  and  the  are  the  values  of 
F  at  the  grid  points  specihed  by  Pierce’s  algorithm.  The  grid  points  are  on 
concentric  circles  in  the  annulus,  with  the  i  index  identifying  the  circles  and 
the  j  index  identifying  the  points  on  a  particular  circle.  Each  circle  has  the 
same  number  of  grid  points,  Auc,  where  Uc  =  the  number  of  circles,  so  the  total 
number  of  grid  points  is  4(nc)^.  Computations  necessary  for  determining  the 
point  locations  and  the  weights  were  done  by  Lowan,  Davids,  and  Levenson 
[24] .  The  result  is  that 

(/)  =  ^  -  tan^  6^-^)  E  E  (12) 

i  j 

This  expression  was  used  to  compute  and  with  equations 

(7a)  and  (8a)  used  to  compute  jstokes  points.  The 

number  of  grid  points  required  for  convergence  becomes  larger  if  the  scat¬ 
tering  particle’s  diameter  is  increased,  because  the  distribution  of  scattered 
light  becomes  more  complex  for  larger  particles. 

3.2  Unpolarized  Incident  Wave 

To  obtain  the  integrands  for  an  unpolarized  incident  wave  most  easily,  we 
average  equations  (7a)  and  (8a)  with  respect  to  0  over  the  interval  0  <  0  < 
27r,  producing 

(V^tokesy  ^  (13a) 

(^I^^oky^  =  Sn.  (13b) 

These  integrands  can  be  used  in  equation  (11)  to  obtain  l^{yStokes'\U\^ 
l^^j-stokes^u Y  Equation  (7b)  shows  that  5'4i  =  0  if  S'3  =  54  =  0.  This  means 
that  there  is  no  false  signal  for  any  aperture  position  when  the  incident  wave 
is  unpolarized. 
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4.  Analytical  Integration  of  Stokes  Parameters 


4.1  Linearly  Polarized  Incident  Wave 

The  initial  step  in  developing  the  analytical  theory  to  be  given  here  is  to 
rewrite  equations  (4)  in  vector  form  as 

yStokes  ^  ^  E*)/{E,nc)^  (14a) 

jStokes  ^  f.2^2^-^  .  (14b) 


Equation  (5)  can  be  substituted  directly  into  equations  (14)  if  scattering 
into  an  aperture  centered  on  the  beam  axis  is  the  only  aperture  location  of 
interest.  Since  we  are  interested  in  computing  and 

desired  aperture  location,  it  is  necessary  to  obtain  expansions  for  Mo-i„  and 
No-in  as  a  linear  combination  of  vector  spherical  harmonics  in  the  aperture 
system.  Such  an  expansion  has  previously  been  given  [16]  for  the  vector 
spherical  harmonic  used  by  Jackson  [20],  and  this  expansion  is  used  here  to 
obtain  expansions  for  M^-in  and  N^-i^. 


The  theory  given  here  is  simplihed  by  the  introduction  of  complex  vector 


spherical  harmonics 


where 


and 


=  Mil  +  iM®„,  (15a) 

flSbi 

^^ran  —  ^^eran  '  ) 


mIL 

—  A 

(16a) 

amn 

=  A  N® 

(16b) 

(2n  +  l)(n  —  m)l 

(16c) 

■mn  —  \ 

2(n  +  m)l 
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is  the  normalizing  factor  for  Legendre  functions.  The  vector  spherical  har¬ 
monics  of  the  beam  system  can  now  be  written  as  a  linear  combination  of 
the  complex  vector  spherical  harmonics  of  the  aperture  system.  A  derivation 
of  this  expansion  for  Mo-in  is  outlined  in  appendix  B,  and  the  expansion  of 
M  a-in  rnay  be  written  as 


m2>„=  (17a) 

m=—n 

where  a  superscript  A  on  vector  spherical  harmonics  indicates  the  use  of  the 
aperture  system  variables  and  basis  vectors  (see  sect.  2.2).  Operating  on  both 
sides  of  equation  (16a)  with  ((1/A;)  v  x)  then  gives 


'  aln 

where  the  matrix  elements  femn  and  fomn  are  dehned  by  the  matrix  equation 


NL1=  E  (17b) 

m=—n 


fe 

fo 


(-l)™+iexp(im7) 


(2n  +  1) 

cos  P) 


I  sm  CK  \  _  ,  /  cos  a 

TT^mtCOSPjl  +r„m(cOS/l) 

—i  cos  a  J  \  sin  a 


(17c) 


To  avoid  an  unnecessary  digression,  we  dehne  and  discuss  the  functions  Unm 
and  Tnm  in  appendix  A. 

Substituting  the  expansions  of  equations  (17)  into  equation  (5)  allows  Bohren’s 
result  to  be  rewritten  as 


cxD  n 


n=l  m=—n 


2n  d- 1 
n{n  +  1) 


rmn  ran  '  rmn  ran 


where 


t/mn  —  ^  emn  dnfomn)  ) 

'^mn  ~  ^nferan  bn  fomn- 


(18a) 


(18b) 

(18c) 
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The  aperture  is  located  in  the  far  held  of  the  sphere.  In  the  far-held  approx¬ 
imation,  Zn  =  h^\kr)  — >  (— gives  the  result  that 


mn 

mn 


„exp(i/cr)  exp(im0"^) 


TTumicOsO^)  \  rnm{cOs9^)  \ 

.  \TnmicOs6^)J  \  Tf^m  (cOS  6'^)  / 


(19) 


Substituting  equations  (18)  into  equations  (14)  and  the  results  into  equation 
(2a)  where  the  integrations  over  0  <  0"^  <  27r  are  performed  gives  the  result 
that 


j^Stokes\^ 

j^jStokes\^ 


^  “  (2n  +  l)(2n'  +  l) 

All  n{n  +  l)n'(n'  +  1) 

(  \  I  (  \ 

1  ^  j.  I  '^nn'm  '  I  I  ^nn'm 

\  ^  .dot  I  \  ^  .cross  I 

^  ^  'rnn'm  '  ^  'rnn'm  '  -• 


(20a) 


where  Af2  =  27r(cos6*An  —  cos6';^j^x)) 


i^nn'm  =  ’  (20^) 

i^n:'Z  ^  Cn^l'  +  ^LCn'  >  (20c) 

=  /  d^FSM.  (21a) 

(21b) 

A^max 

where  /Xmin  =  cos{6^^J,  /Xmax  =  cos(6'i^^^),  and 

■^nn^m  (/^)  ~  “1“  ^nm(/^)^n^m(/^)*  (^1^) 
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^aoz  ^aoz  ■.  •  /nr\  \  1  'J  i 

Since  and  equation  (20a)  can  be  rewritten 

as 


I^Stokes\^ 

^jStokes^ 


71 


oo  oo  min(n,n') 

EE  E 


n=l  n'=l  m=0 


Nm{2n  +  l){2n'  +  1) 
n{n  +  l)n'{n'  +  1) 


{+)cross 

nn'm 

{-\-)dot 

nn'm 


jdot 

^nn'm 


+ 


^nn'm  \  jCross 
{-)cross  I  nn'm 


(22a) 


where  N^,  is  the  Neumann  factor 


Nm  =  1  if  m  =  0, 

Nm  =  2  if  m  Y  0; 

(22b) 

^(±)dot  _  I  dot  1  idot 

^  nn'm  —  rnn'm  rnn'^—m’’ 

(22c) 

^{±)cross  _  1  cross  i  i  cross 
^  nn'm  —  rnn'm  ^  rnn'^—m' 

(22d) 

Equations  (22)  can  be  programmed  to  compute  the  most  general  case  for 
scattering  of  a  linearly  polarized  incident  wave. 


4.2  Unpolarized  Incident  Wave 

We  obtained  the  solutions  for  an  unpolarized  incident  wave  from  equations 
(22)  by  averaging  over  angle  a  to  produce 


Stokes 


U 


^  ^  ^  minin')  iv^(2n  +  1) (2n^  +  1) (-1) 

^0  n{n  +  l)n'(n'  +  1) 

I  (^(-)dot\^  ^cross' 

\^nn'm  I  ^nn'm  nn'm  I  ^nn'm  ’ 


(23a) 


jStokes\  \  ^ 


77  minp^n')  iV^(2n  +  1)  (2n^  +  1) 

n=in'=i  ho  n(n  +  l)n'(n'  +  l) 


( ^ -wdot 
nn'm  J  ^nn't 


+ 


nn'r, 


) 


U-c 

i! 


(23b) 
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where 


2tt 


Averaging  the  functions  with  respect  to  angle  a  gives 

^nn'm(cOS/3), 


f^(+)cr-oss\^  _ 

4 

1  r 
z- 

\  nn  m  J 

_(2n  +  l)(2n'  +  1)_ 

2  L 

^(+)cross  ^(-\-)cross* 

^nn'  ^n'n 


4 

.1  r 

i  — 

\  nn  m  J 

_(2n  +  l)(2n'  +  1)_ 

2  L 

\  nn 


{+)dot\  ^ 


^{-)dot  _  ^{-)dot*^  uncross  f  n\ 

^nn  ^n'n  ^  nn' P  )  i 


(2n  +  l)(2n'  +  1) 


(^(  Y 

\  nn  r 


U 


where 


(2n  +  l)(2n'  +  1) 


(C«--^  +  Ci7i™)F:X(cos/3), 


Citr^°^^  =  2{ar.c:,  +  d^bY), 

dJ^^^2{dna:,  +  bncY), 

=  (flnfln'  +  bnbY  +  C„C*/  +  dndY)  , 

ci-yross  ^  2  (a  J>l,-dncY)- 


(23c) 


(24a) 

(24b) 

(24c) 

(24d) 

(25a) 

(25b) 

(25c) 

(25d) 


Interchanging  n  and  n'  (a  standard  symmetry  operation)  in  the  C'*/^  terms 
and  noting  that  F^X(cos/3),  F^°)^(cos/3),  and  I„“,^  are  unchanged 

after  this  operation  allows  equations  (23)  to  be  rewritten  as 
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u 


Stokes 


Att 


An 


oo  oo  min(n,n') 

EE  E 

n=l  n'=l  m=0 


X 


N„  [lm(cW°’'°")  F„a.(cos  /3)C„.  +  Im(ci7/°‘)  FS‘l(a»  /3)C" 

n{n  +  l)n'(n'  +  1) 


^cros& 
rr 


jStokes^^ ^ 


oo  oo  min(n,nO 


AO 

n=ln'=l  m=0 


n{n  +  l)n'(n'  +  1) 


,  (26a) 


X 


.  (26b) 


5.  Reduction  for  Scattering  into  a  Solid  Angle 
47r  Steradians 


We  have  investigated  experimental  configurations  collecting  almost  all  the 
scattered  light.  The  light  is  scattered  into  a  solid  angle  Ahl  of  approximately 
dTT  sr  (i.e.,  =  0  and  9^^^  =  n  radians  =  180°).  The  average  Stokes 

parameters  can  be  calculated  for  this  case  without  the  added  complication 
of  a  beam  coordinate  system  and  Euler  angles.  However,  the  mathematical 
technique  required  for  the  simplification  of  equations  (22)  and  (26)  when 
AH  =  An  sr  may  be  of  interest  to  theorists  who  use  angular  momentum 
theory.  The  result  must  be  independent  of  the  Euler  angles  since  all  the 
scattered  light  is  collected.  Further,  the  results  for  linearly  polarized  and 
unpolarized  incident  waves  must  be  identical. 


5.1  Linearly  Polarized  Incident  Wave 

Mathematically,  the  simplification  of  equations  (22)  and  (26)  occurs  because 

'^dot  ■|A‘min  =  l  i  \  r  1  1  r\  •  T  j  1  1  An 

Ln'm  =  n[n+  IjOnn'  and  i-rin'm  =  U  as  indicated  by  Artken 

f^max—  1  ^  A^max—  1 

[22]  and  appendix  C.  The  result  is  that 


=  V 


{2n+iy 


4n(n  +  1) 


I  -‘-y  rn^  nnrr. 


oo 

/  jStokes\  _  \  ^ 

A  ^\An=4n  ~  ^ 


Tl  4n(n  +  1)  ' 


(27b) 


With  equation  (D-3)  derived  in  appendix  D,  it  is  not  difficult  to  show  that 

(-1)  Y.  lo„<  +  ,  (28a) 

i:  =  1*1  [l“«t  +  linf"  +  \Cnf  +  .  (28b) 

m=0  ^  “I"  -‘-1 
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These  expressions  give 


An.4.  =  +  !)■■“  I“”<  +  • 


n=l 


^Stokes'^  _  -  ^(2n+  1)  \anf  +  \bnf  +  \Cnf  +  \d„ 

^  n=l 


(29a) 

(29b) 


These  are  the  average  Stokes  parameters  for  scattering  from  a  linearly  po¬ 
larized  incident  wave  into  dvr  sr. 


The  scattering  efficiency  (^scattering  is  connected  to  (^i Stokes'^ 
see  by  comparing  their  dehnitions  and  observing  that 

4 

_  I  ^  j  sjLtjf\ic.a 

'  J  AO=47r 


,  as  we  can 


^  _  _  /  jbtokes\ 

scdtteving  /t  no  \  / 


This  gives  the  result  that 

2  “ 

(ka) 


Q  scattering  {u^\2 


{ka) 


^(2n-|-l)  \an\  +  |&n|  +  |Cn|  +  |d„ 


(30a) 


(30b) 


If  there  is  no  optical  activity,  then  a„  =  bn  =  and  Cn  =  dn  =  0,  so 
that  the  conventional  expression  for  Mie  scattering  efficiency  is  obtained. 


Unpolarized  Incident  Wave 

For  the  unpolarized  wave  scattering  into  Af2  =  47r  sr, 


jdot 
^nn'm 

^cros 


■^cross 

nn'm 


A^min  —  1 
A^max—  1 
A^min  —  1 
A^max—  1 


n{n  +  l)Snn' 
0 


and 


give 


Stokes 


U 


AQ=47r 


+  1) 


m=0 


jStokes\^ 


U 


An=47r 


oo  "Rp  ( 

V  y^nn  J 

n(n  +  1) 


m=0 


(31a) 

(31b) 


With  the  expressions  given  in  appendix  D,  it  is  found  that 


so  that 


n  ^ 

^mi^nnm(cos/3)  =  -n{n  +  l)(2n  +  1), 

m=0  ^ 

1 1  1 

^  =-i:(2n  +  l)  Im(cg>™>), 

An  4-it  Z 

T I  1 

_  =-i:(2n  +  l)Re(cW“), 

An  4ir  n=l 


showing  that 


(32) 

(33a) 

(33b) 

(34a) 

(34b) 


19 


6.  False  Signal 


In  sections  3.1  and  3.2,  we  show  that  the  term  “false  signal”  (of  optical 
activity)  is  eqnivalent  to  a  nonzero  measnrement  of  for  scattering 

from  an  optically  inactive  particle  for  which  c„  =  —dn  =  0. 

For  a  linearly  polarized  incident  wave  (sect.  3.1),  we  show  that  there  is  no 
false  signal  for  aperture  locations  such  that  a  =  0  or  7r/2;  i.e.,  =  Q 

for  scattering  from  a  spherical  particle  with  no  optical  activity  when  the 
aperture  is  centered  in  the  x-z  plane  or  in  the  y-z  plane.  We  can  also  obtain 
this  result  by  analyzing  equation  (22a).  If  a  =  0  or  7r/2,  or  if  /?  =  0  or  n, 
then  =  0  optically  inactive  spherical  particles,  so  that 

equation  (22a)  gives  the  result  that  =  Q. 

For  an  unpolarized  incident  beam,  =  0  for  optically  inac¬ 

tive  spheres,  so  that  equation  (26a)  gives  the  result 

section  3.2,  our  results  show  that  there  is  no  false  signal  for  an  unpolarized 
incident  beam. 
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7.  Numerical  Calculations 


In  this  section,  we  give  the  results  of  some  computations  for  the  beneht  of 
computational  researchers  who  wish  to  verify  computer  programs  that  they 
develop  with  the  theory  given  here.  Using  Mie  theory,  Bohren  and  Huffman 
give  computed  values  of  and  6^*®  for  Re(iV)  =  1.33,  Im(iV)  =  1.0(10“®), 
and  the  Mie  size  parameter  X  =  ttD/X  =  3.0.  All  the  computations  given 
here  are  for  an  optically  active  sphere  with  X  =  3.0,  Re(A’i)  =  1.33, 
Im(A'i)  =  2.0(10“®),  Re{Nji)  =  1.33,  and  =  0.  For  these  values, 

the  resulting  computed  values  of  a„,  bn,  and  c„  are  given  in  table  1. 

7.1  Aperture  Centered  on  the  Beam  Axis 

In  tables  2,  3,  and  4,  computations  of  ^jStokes'^ 

given  for  scattering  into  a  sensor  aperture  centered  on  the  beam  axis  {a  = 
(3  =  0°).  For  this  particular  aperture  location,  the  results  are  identical  for 
polarized  and  unpolarized  incident  beams. 

Table  2  is  for  an  aperture  with  beam  block  for  which  dmm  =  5°  and  0max  = 
60°.  We  did  these  computations  with  the  numerical  integration  scheme,  giving 


Table  1.  Computed  values  of  a„,  6„,  and  c„. 


n 

R’^(O'n) 

Im(a„) 

Re(6„) 

Im(&„) 

1 

0.516306(10°) 

-0.499734(10°) 

0.737672(10°) 

-0.439900(10°) 

2 

0.341921(10°) 

-0.474353(10°) 

0.400793(10°) 

-0.490059(10°) 

3 

0.484668(10“!) 

-0.214750(10°) 

0.935528(10“2) 

-0.962692(10“!) 

4 

0.103458(10“2) 

-0.321483(10“!) 

0.688105(10“'!) 

-0.829491(10“°) 

5 

0.903755(10“°) 

-0.300622(10“2) 

0.283089(10“°) 

-0.532042(10“®) 

n 

Im(c„) 

1 

-0.265630(10“!") 

-0.717984(10“®) 

2 

-0.261447(10“!") 

0.697491(10“®) 

3 

-0.202011(10“®) 

0.611957(10“®) 

4 

-0.268083(10“!°) 

0.662132(10“°) 

5 

-0.186204(10“!2) 

0.526254(10“!°) 
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Table  2.  ,  ^yStokes^^  ^Stokes^  for  a  =  /3  =  0,  6»min  =  5°,  and  e^ax  =  60°; 

computed  with  numerical  integration  method  with  identical  results  for  polarized  and  un¬ 
polarized  incident  waves. 


Number 
of  circles 

Number 
of  points 

peffective 

I^Stokes\^ 

^jStokes'j 

2 

16 

0.534274(10-'^) 

0.540146(10-®) 

0.101099(102) 

4 

64 

0.499270(10-'^) 

0.700091(10-®) 

0.140223(102) 

6 

128 

0.496561(10-'^) 

0.708648(10-®) 

0.142711(102) 

8 

256 

0.496438(10-'^) 

0.708921(10-®) 

0.142802(102) 

10 

400 

0.496434(10-'^) 

0.708928(10-®) 

0.142804(102) 

12 

576 

0.496434(10-'^) 

0.709932(10-®) 

0.142805(102) 

14 

784 

0.496434(10-'^) 

0.708928(10-®) 

0.142804(102) 

Table  3.  ^yStokes^^  ^Stokes^  for  a  =  P  =  0,  =  0,  and  6»„,ax  as  indicated; 

computed  with  analytical  method  with  identical  results  for  polarized  and  unpolarized 
incident  waves. 


Ora.,  n 

peffective 

l^jStokes'j 

60 

0.495059(10-'^) 

0.717150(10-®) 

0.144861(102) 

120 

0.516809(10-'^) 

0.268489(10-®) 

0.519513(101) 

180 

0.517047(10-'^) 

0.203983(10-®) 

0.394514(101) 

results  for  an  increasing  number  of  points  to  demonstrate  the  convergence  of 
this  method.  The  results  given  in  the  bottom  row  of  table  2  were  found  to 
agree  exactly  with  those  given  by  the  analytical  method. 

In  table  3,  computed  values  of  p^ffecuve^  (y Stokes'^  ^  i^jstokes\^  given  for 

6'min  =  0°  and  6*max  as  indicated.  These  computations  were  made  with  the 
analytical  method.  The  values  of  table  3  for  =  60°  were  verified  with 
the  numerical  integration  method,  and  the  computations  for  6'inax  =  180° 
were  verihed  with  the  specialized  analytical  result  for  scattering  into  47r  sr. 
Exact  agreement  was  obtained  in  each  case. 

Off-Axis  Aperture  Locations 

An  off-axis  aperture  location  was  chosen  with  a  =  30°  and  (3  =  60°.  In  table  4, 
computations  of  P^ffecUve^  yStokes\^^  l^jStokes\^  given  for  6'min  =  0°  and 

^max  as  indicated.  The  computations  of  table  4a  are  for  a  polarized  incident 
wave,  and  the  computations  of  table  4b  are  for  an  unpolarized  incident  wave. 


Table  4.  (y stokes-^ ^  ^Stokes'^  £qj,  incident  wave  scattering  into  aperture 

with  a  =  30°,  (3  =  60°,  0min  =  0°,  and  ^max  as  indicated;  computed  with  analytical 
method:  (a)  polarized  and  (b)  unpolarized. 


(a) 


0max  (°) 

p  effective 

^'Y'Stokes'^ 

^jStokes'j 

60 

-0.434718(10-2) 

-0.298083(10-1) 

0.685693(101) 

120 

-0.223394(10-2) 

-0.114976(10-1) 

0.514676(101) 

180 

0.517047(10-'^) 

0.203983(10-6) 

0.394514(101) 

(b) 

0max  (°) 

p  effective 

^'Y'Stokes^ 

^jStokes'j 

60 

0.508052(10-'^) 

0.349986(10-6) 

0.688877(101) 

120 

0.512351(10-'^) 

0.263510(10-6) 

0.514315(101) 

180 

0.517047(10-'^) 

0.203983(10-6) 

0.394514(101) 

The  values  in  these  tables  were  computed  with  a  Cray  and  a  desktop  per¬ 
sonal  computer.  All  results  given  by  these  two  computers  agreed  except  for 
the  values  of  p^ifecawe  iyStokes\^  obtained  for  table  4a  for  6*jnax  =  180°.  The 
results  given  by  the  personal  computer  were  incorrect  in  this  case  because 
of  a  lack  of  precision.  Since  many  computations  are  performed  on  personal 
computers,  we  considered  it  important  to  explain  this  discrepancy.  The  spe¬ 
cialized  result  for  obtained  for  scattering  into  47r  sr  follows  from  the 

observation  that 

in  n  n 

--  ^  =  E  iV^|/e„^npIm(a„c;)  +  E  iV™|/on^npIm(d„6:) 

m=0  m=0  m=0 

n 

-  E  N^^e[femnf*omn][^^i.dnCl)  +  Im(an&;)].  (35) 

m=0 


The  difficulty  arises  when  Yl'in=o  -^mRe  [/emn/omn]  exactly  zero  because 

this  factor  is  multiplied  by  Im(a„6* )  where  a„  and  bn  are  much  larger  than 
dn  =  —Cn-  Although  we  found  theoretically  that  A^Re  [/emn/omn]  = 

0,  the  limited  precision  of  personal  computers  is  such  that  if  each  term  of 
Re  [/emn/omn]  not  exactly  zero,  then  the  sum  differs  from  zero  enough  to 
give  an  erroneous  result  when  multiplied  by  Im(a„6*).  The  difficulty  does  not 
arise  when  high-precision  machines  like  Cray  computers  are  used  or  when  a 
personal  computer  is  used  with  the  factor  Re  [/emn/omn]  =  0  in  each  term. 
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This  factor  is  zero  if  a  =  0  or  7r/2,  and  it  is  also  zero  if  /3  =  0  or  tt.  Using  these 
aperture  positions  corresponds  to  elimination  of  the  false  signal  as  discussed 
in  section  6. 
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8.  Summary 


The  objective  of  the  research  initiated  earlier  [11]  was  to  develop  a  method  for 
identifying  optically  active  spheres  by  analyzing  light  scattered  from  a  laser 
beam  into  a  circular  sensor  aperture.  We  claimed  that  this  could  be  done  by 
computing  the  effective  degree  of  circular  polarization  for  apertures  located 
in  certain  positions  that  generate  no  false  signal  of  optical  activity,  and  these 
positions  were  identihed.  Here,  we  prove  this  claim,  and  two  methods  for 
performing  the  necessary  computations  are  developed  and  discussed. 

We  used  Bohren’s  theory  for  light  scattering  by  an  optically  active  sphere  as 
a  starting  point  to  develop  the  two  methods  given  here,  which  average  Stokes 
parameters  over  a  circular  sensor  aperture  located  in  any  desired  direction 
from  the  scattering  sphere.  The  effective  degree  of  circular  polarization  is 
dehned  as  the  ratio  of  two  averaged  Stokes  parameters. 

One  method  is  based  on  a  two-dimensional  numerical  integration  scheme  that 
explicitly  shows  the  connection  with  Mueller  matrix  elements;  this  method  is 
referred  to  as  the  numerical  integration  method.  The  other  method  is  based 
on  expansions  of  vector  spherical  harmonics  and  is  referred  to  as  the  analyti¬ 
cal  method.  The  advantages  and  disadvantages  of  each  method  are  discussed. 
For  a  dedicated  programmer  with  a  fast  computer,  the  analytical  method  is 
preferred  since  it  can  treat  apertures  of  any  size  without  convergence  checks. 
Either  method  can  treat  both  polarized  and  unpolarized  incident  beams. 
Equations  are  derived  that  allow  reduction  of  the  general  analytical  expres¬ 
sions  in  the  special  case  of  scattering  into  dvr  sr.  Several  tables  of  numerical 
results  are  given  to  facilitate  the  verihcation  of  computer  codes,  and  we  dis¬ 
cuss  potential  computational  difficulties  with  personal  computers  because  of 
a  lack  of  precision.  We  pointed  out  that  these  difficulties  do  not  arise  when 
the  aperture  is  located  in  a  position  that  generates  no  false  signal.  Mathemat¬ 
ical  details  related  to  special  functions,  expansions,  integrals,  and  reductions 
are  given  in  four  appendices. 
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Appendix  A.  Definition  and  Computation  of  the  iinm 
and  Tnm  Functions 


The  TTnmil^)  and  Tnmil^)  functions  are  defined  here  as 

TTnmlh)  =  ,  (A-la) 

vi  -/^ 

^nm(h')  =  h  (*^/ (A-lb) 

where  Pnm.{.lj)  is  the  normalized  associated  Legendre  function  as  dehned  by 
Belousov  [25]  so  that 


(A-lc) 


where  Pnm{lA  =  (1  ~  is  fhe  unnormalized  associated  Legen¬ 

dre  function  as  dehned  by  Arfken  [22]  and  Stratton  [26]  and  used  in  standard 
light  scattering  references  by  Bohren  and  Huffman  [1]  and  van  de  Hulst  [23]. 
These  Legendre  functions  differ  by  a  factor  of  (—I)’”  from  the  Legendre  func¬ 
tions  used  by  Jackson  [20],  Abramowitz  and  Stegun  [27],  and  Gradshteyn  and 
Ryzhik  [28].  These  last  two  references  are  of  particular  interest  for  their  ex¬ 
tensive  recursion  relations,  but  care  must  be  exercised  because  of  the  (— I)™' 
factor  difference.  The  Unm  and  Tnm  functions  are  dehned  here  to  agree  with 
Fuller’s  dehnitions  [29]  except  with  an  additional  factor  of  Amn.  There  is 
a  computational  advantage  to  using  normalized  Legendre  functions  because 
the  factorials  in  A^n  do  not  have  to  be  computed  separately  but  are  easily  in¬ 
corporated  into  the  recursion  relations.  Also,  for  negative  values  of  m,  the  re¬ 
sulting  functions  are  Pn-m(h)  =  (-l)™Pnm(h),  -m(h)  =  (“ (h) , 

and  Tn-mi/J')  =  (— These  relations  are  simpler  than  those  using 
the  unnormalized  Legendre  functions. 

We  can  conveniently  conpute  the  Tnm(h)  und  functions  using  relations 

obtained  by  rewriting  equation  (8.733.1)  given  by  Gradshteyn  and  Ryzhik 
[28]  as 
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{n  -  m  +  l){n  +  m)Pn,m-i{fj) 


{±}  J{n  +  m  +  l){n-  m)Pn,m+i{l^) 


If  /i  =  1,  then  P„m(l)  =  y (2n  +  l)/25m,o  gives  the  result 


7rnm(l)  1  _  1  /(2n  + l)n(n  + 1) 


{±}  Sm-l)- 


rnm(l)  J  2  V  2 

It  is  helpful  to  realize  that  this  expression  may  be  rewritten  as 


^  nm  ( 1 ) 
V  nm ( 1 ) 


Ai^nnjn  +  1) 

2 

_  (2n  +  1) 
4^1, n 


{Sm,l  {±}  Sm-l) 


{6m, 1  {±}(5  171,-1)* 


(A-2a) 


(A-2b) 


(A-2c) 


Appendix  B.  Expansion  of  the  Vector  Spherical 
Harmonics 

The  expansion  of  Jackson’s  vector  spherical  harmonic  can  be  written  as  [16] 

n 

x„„=  (B-la) 

m=—n 

where  0),  X;^^  =  0^),  and 

=  exp{ima)dM^^{p)  exp(im7).  (B-lb) 


_ j - 

Comparing  the  definitions  of  and  X„m  gives  =  J2nn{n  +  1) 

{—l)^~^^iZnXnm  SO  that  eqnation  (B-la)  can  be  used  to  obtain  the  relation 


mS,„=  E  Dri,™(a/37)(-l)”+‘M27 

m=—n 

(B-2a) 

where 

=  exp(±ia)d[^i  ,„(0)  exp(im7), 

(B-2b) 

=  \  f  ,  TVO  ,  a  jTTnm  (cos  0)  (cos  0)]. 

V  n{n  +  l)(2n  J-  1) 

(B-2c) 

The  equations 

=  [m®  +  (-1)”^M^^^  1  /2, 

(B-3a) 

-  (-1)”^M^^^  1  /(20 

(B-3b) 

(with  m  =  1),  and  equation  (B-2)  then  give 


n 

=  S  UXn,  (B-4) 

m=—n 

which  can  be  rewritten  as  equations  (17a)  and  (17b)  with  the  expansion 
coefficients  famn  given  by  equation  (17c). 
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Appendix  C.  Evaluation  of  Ijl'm  K 


The  theory  of  evaluating  the  integrals  Inn'm  Inn'm  mainly  done  by 
Chu  and  Robinson  [15]  with  a  later  contribution  by  Son,  Farmer,  and  Giel 
[18].  The  results  are  listed  here  in  the  notation  of  this  report  for  the  conve¬ 
nience  of  the  reader,  with  a  few  comments  concerning  proof  of  the  results. 

By  differentiating  both  sides  of  the  equation,  we  easily  show  that 

CC  =  [-mPnMPnM]  7"  .  (C-1) 

For  what  follows,  it  is  helpful  to  rewrite  the  differential  equation  for  the 
associated  Legendre  function  as 

^  (^Y^l  -  =  n{n  +  l)Pnm(/i)  “  '  (^"2) 

With  this  result,  we  can  easily  show  (by  differentiating  both  sides)  that  if 
n  ^  n',  then 


Equation  (C-4c)  can  be  proved  with  the  help  of  a  relation,  = 

(/i7f„m(h)  +  /  J {n  +  m){n  —  m  +  1) ,  obtained  from  equation  (A-2a). 
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Appendix  D.  Equations  Giving  Reduction  for 
Scattering  into  An  Steradians 


Equations  required  to  demonstrate  the  simplification  that  occurs  for  scat¬ 
tering  into  dTT  sr  can  be  obtained  in  the  following  way:  if  a  =  7  =  0  and 
^  0  and  if  6*  =  0  =  0,  then  6^  =  (3  and  (j)^  =  tt  while  6  =  —6  =  x  and 
cj)  =  —0  =  y.  Substituting  this  set  of  conditions  into  equation  (B-4)  and 

using  =  Ainn{n  -|-  l)/2  from  equation  (A-2c)  gives 

n 

XI  =  X  [^nm(cos/3)]^  = -n(n-h  l)(2n-h  1),  (D-la) 

m=—n  rn=—n 

n 

X/  '^nm  (cos /3)t77772  (cos /3)  —  0.  (D-lb^ 


Since  and  Tn-m{y)  =  nm{y) ,  h  is  easy  to 

show  that  fa-m,n  =  famn'  With  this  expression  and  equation  (D-3),  it 

is  easy  to  obtain  equations  (28). 


31 


References 


1.  C.  F.  Bohren  and  D.  R.  Huffman,  Absorption  and  Scattering  of  Light 
by  Small  Particles  (Wiley-Interscience,  New  York,  1983),  pp.  45,  49,  63, 
185-190. 

2.  G.  Holzwarth,  D.  J.  Gordon,  J.  E.  McGinnes,  B.  P.  Dorman,  and  M.  F. 
Maestre,  “Mie  scattering  contributions  to  the  optical  density  and  circu¬ 
lar  dicroism  of  t2  bacteriophage,”  Biochemistry  13,  126-131  (1974). 

3.  D.  J.  Gordon,  “Mie  scattering  by  optically  active  particles,”  Biochem¬ 
istry  11,  413-420  (1972). 

4.  K.  D.  Lofftus,  M.  S.  Quinby-Hunt,  A.  J.  Hunt,  F.  Livolant,  and  M. 
Maestre,  “Light  scattering  by  Prorocentrum  micans:  a  new  method  and 
results,”  Appl.  Opt.  31,  2924-2931  (1992). 

5.  B.  V.  Bronk,  S.  D.  Druger,  J.  Gzege,  and  W.  P.  van  de  Merwe,  “Mea¬ 
suring  diameters  of  rod-shaped  bacteria  in  vivo  with  polarized  light 
scattering,”  Biophys.  J.  69,  1170-1177  (1995). 

6.  G.  F.  Bohren,  “Light  scattering  by  an  optically  active  sphere,”  Ghem. 
Phys.  Lett.  29,  458-462  (1974). 

7.  G.  F.  Bohren,  “Scattering  of  electromagnetic  waves  by  an  optically  ac¬ 
tive  spherical  shell,”  J.  Ghem.  Phys.  62,  1566-1571  (1975). 

8.  G.  F.  Bohren,  “Scattering  of  an  optically  active  cylinder,”  J.  Golloid 
Interface  Sci.  66,  105-109  (1978). 

9.  M.  S.  Kluskens  and  E.  H.  Newman,  “Scattering  by  a  multilayer  chiral 
cylinder,”  IEEE  Trans.  Antennas  Propag.  39,  91-96  (1991). 

10.  M.F.R.  Gooray  and  1.  R.  Giric,  “Wave  scattering  by  a  chiral  spheroid,” 
J.  Opt.  Soc.  Am.  A  10,  1197-1203  (1993). 

11.  D.  L.  Rosen  and  J.  D.  Pendleton,  “Detection  of  biological  particles 
by  the  use  of  circular  dichroism  measurements  improved  by  scattering 
theory,”  Appl.  Opt.  34,  5875-5884  (1995). 


33 


12.  W.  H.  Pierce,  “Numerical  integration  over  the  planar  annulus,”  J.  Soc. 
Ind.  Appl.  Math.  5,  66-73  (1957). 

13.  P.  Chylek,  “Mie  scattering  into  the  backward  hemisphere,”  J.  Opt.  Soc. 
Am.  63,  1467-1471  (1973). 

14.  W.  J.  Wiscombe  and  P.  Chylek,  “Mie  scattering  between  any  two  an¬ 
gles,”  J.  Opt.  Soc.  Am.  67,  572-573  (1977). 

15.  W.  P.  Chu  and  D.  M.  Robinson,  “Scattering  from  a  moving  spherical 
particle  by  two  crossed  coherent  plane  waves,”  Appl.  Opt.  16,  619-626 
(1977). 

16.  J.  D.  Pendleton,  “Mie  scattering  into  apertures,”  J.  Opt.  Soc.  Am.  72, 
1029-1033  (1982). 

17.  J.  D.  Pendleton,  A  Generalized  Mie  Theory  Solution  and  Its  Application 
to  Particle  Sizing  Interferometry,  PhD  thesis  (University  of  Tennessee, 
Knoxville,  1982),  p.  90. 

18.  J.  Y.  Son,  W.  M.  Farmer,  and  T.  V.  Giel,  Jr.,  “New  optical  geometry  for 
the  particle  sizing  interferometer,”  Appl.  Opt.  25,  4332-4337  (1986). 

19.  J.  Y.  Son,  Multiple  Methods  for  Obtaining  Particle  Size  Distribution 
With  a  Particle  Sizing  Interferometer,  PhD  thesis  (University  of  Ten¬ 
nessee,  Knoxville,  1985). 

20.  J.  D.  Jackson,  Classical  Electrodynamics,  2nd  ed.  (Wiley,  New  York, 
1975),  p.  746. 

21.  M.  E.  Rose,  Elementary  Theory  of  Angular  Momentum  (Wiley,  New 
York,  1957),  Chap.  4. 

22.  G.  Arfken,  Mathematical  Methods  for  Physicists,  3rd  ed.  (Academic 
Press,  San  Diego,  1985),  pp.  198-200,  253,  678. 

23.  H.  C.  van  de  Hulst,  Light  Scattering  by  Small  Particles  (Dover,  New 
York,  1981),  pp.  34,  124. 

24.  A.  N.  Lowan,  N.  Davids,  and  A.  Levenson,  “Table  of  the  zeros  of  the  Leg¬ 
endre  polynomials  of  order  1-16  and  the  weight  coefficients  for  Gauss’ 
mechanical  quadrature  formula,”  Am.  Math.  Soc.  Bull.  48,  739-743 
(1942);  “Errata,”  49,  939  (1943). 


34 


25.  S.  L.  Belousov,  Tables  of  Normalized  Associated  Legendre  Polynomials 
(Pergamon,  New  York,  1962). 

26.  J.  A.  Stratton,  Electromagnetic  Theory  (McGraw-Hill,  New  York,  1941), 
p.  401. 

27.  M.  Abramowitz  and  I.  A.  Stegun,  Handbook  of  Mathematical  Functions 
(Dover,  New  York,  1964). 

28.  I.  S.  Gradsliteyn  and  I.  M.  Ryzhik,  Table  of  Integrals,  Series,  and  Prod¬ 
ucts  (Academic,  New  York,  1980),  pp.  1005,  1008. 

29.  K.  A.  Fuller,  “Scattering  and  absorption  cross  sections  of  compounded 
spheres.  I.  Theory  for  external  aggregation,”  J.  Opt.  Soc.  Am.  A  11, 
3251-3260  (1994). 


35 


Distribution 


Admnstr 

Defns  Techl  Info  Ctr 
Attn  DTIC-OCP 

8725  John  J  Kingman  Rd  Ste  0944 
FT  Belvoir  VA  22060-6218 

Central  Intllgnc  Agency  Dir  DB  Standard 
Attn  GE  47  QB 
Washington  DC  20505 

Chairman  Joint  Chiefs  of  Staff 
Attn  J5  R&D  Div 
Washington  DC  20301 

Defns  Intllgnc  Agency 
Attn  DT  2  Wpns  &  Sys  Div 
Washington  DC  20301 

Dir  of  Defns  Rsrch  &  Engrg 
Attn  DD  TWP 
Attn  Engrg 
Washington  DC  20301 

Ofc  of  the  Dir  Rsrch  and  Engrg 
Attn  R  Menz 
Pentagon  Rm  3E1089 
Washington  DC  20301-3080 

Ofc  of  the  Secy  of  Defns 
Attn  ODDRE  (R&AT)  G  Singley 
Attn  ODDRE  (R&AT)  S  Gontarek 
The  Pentagon 

Washington  DC  20301-3080 

OIR  CSB  CRB 
Attn  A  M  Jones 
RB  1413  OHM 
Washington  DC  20505 

OSD 

Attn  OUSD(A&T)/ODDDR&E(R)  J  Lupo 
The  Pentagon 

Washington  DC  20301-7100 

Commanding  Officer 
Attn  NMCB23 
6205  Stuart  Rd  Ste  101 
FT  Belvoir  VA  22060-5275 

CECOM 

Attn  PM  GPS  COL  S  Young 
FT  Monmouth  NJ  07703 


CECOM  RDEC  Elect  System  Div  Dir 
AttnJ  Niemela 
FT  Monmouth  NJ  07703 

CECOM 

Sp  &  Terrestrial  Commctn  Div 
Attn  AMSEL-RD-ST-MC-M  H  Soicher 
FT  Monmouth  NJ  07703-5203 

Dir  of  Chem  &  Nuc  Ops  DA  DCSOPS 
Attn  Techl  Lib 
Washington  DC  20301 

Dpty  Assist  Secy  for  Rsrch  &  Techl 
AttnSARD-TT  F  Milton  Rm  3E479 
The  Pentagon 

Washington  DC  20301-0103 

Hdqtrs  Dept  of  the  Army 
AttnDAMO-FDT  D  Schmidt  Rm  3C514 
400  Army  Pentagon 
Washington  DC  20301-0460 

US  Army  Engrg  Div 
Attn  HNDED  ED 
PO  Box  1500 
Huntsville  AL  35807 

US  Army  ERDEC 
Attn  I  Sindoni 
AttnJ  Embury 
Attn  M  Milham 
Attn  S  Christesen 
Attn  S  Godoff 
Attn  B  Bronk 

Aberdeen  Proving  Ground  MD  21005-5423 

US  Army  Matl  Cmnd 
Dpty  CG  for  RDE  Hdqtrs 
Attn  AMCRD  BG  Beauchamp 
5001  Eisenhower  Ave 
Alexandria  VA  22333-0001 

US  Army  Matl  Cmnd 
Prin  Dpty  for  Acquisition  Hdqrts 
Attn  AMCDCG-A  D  Adams 
5001  Eisenhower  Ave 
Alexandria  VA  22333-0001 


37 


Distribution  (cont'd) 


US  Army  Matl  Cmnd 
Prin  Dpty  for  Techlgy  Hdqrts 
Attn  AMCDCG-T  M  Fisette 
5001  Eisenhower  Ave 
Alexandria  VA  22333-0001 

US  Army  Mis  &  Spc  Intllgnc  Ctr 

Attn  AIAMS  YDL 

Redstone  Arsenal  AL  35898-5500 

US  Army  NGIC 
Attn  Rsrch  &  Data  Branch 
220  7th  Stret  NE 
Charlottesville  VA  22901-5396 

US  Army  Nuc  &  Cheml  Agency 
7150  Heller  Loop  Ste  101 
Springfield  VA  22150-3198 

US  Army  Rsrch  Ofc 
Attn  AMXRO-ICA  B  Mann 
PO  Box  12211 

Research  Triangle  Park  NC  27709-2211 

US  Army  Strtgc  Defns  Cmnd 
Attn  CSSD  H  MPL  Techl  Lib 
Attn  CSSD  H  XM  Davies 
PO  Box  1500 
Huntsville  AL  35807 

US  Military  Academy 
Dept  of  Mathematical  Sci 
Attn  MAJ  D  Engen 
West  Point  NY  10996 

USAASA 

Attn  MOAS-AI  W  Parron 
9325  Gunston  Rd  Ste  N319 
PT  Belvoir  VA  22060-5582 


GPS  Joint  Prog  Ofc  Dir 

Attn  COL  J  Clay 

2435  Vela  Way  Ste  1613 

Los  Angeles  APB  CA  90245-5500 

Special  Assist  to  the  Wing  Cmndr 
Attn  50SW  /  CCX  Capt  P  H  Bernstein 
300  O'Malley  Ave  Ste  20 
Palcon  APB  CO  80912-3020 

DARPA 

Attn  B  Kaspar 

Attn  L  Stotts 

Attn  Techl  Lib 

3701  N  Pairfax  Dr 

Arlington  VA  22203-1714 

US  Dept  of  Energy 
Attn  KK  22  K  Sisson 
Attn  Techl  Lib 
Washington  DC  20585 

ARL  Electromag  Group 

Attn  Campus  Mail  Code  P0250  A  Tucker 

University  of  Texas 

Austin  TX  78712 

US  Army  Rsrch  Lab 

Attn  AMSRL-CI-LL  Techl  Lib  (3  copies) 

Attn  AMSRL-CS-AL-TA  Mail  &  Records 
Mgmt 

Attn  AMSRL-CS-AL-TP  Techl  Pub  (3  copies) 
Attn  AMSRL-IS-EE  D  L  Rosen  (15  copies) 
Attn  AMSRL-IS-EE  J  B  Gillespie  (15  copies) 
Attn  AMSRL-IS-EE  J  D  Pendleton  (15  copies) 
Adelphi  MD  20783-1197 


Chief  of  Nav  OPS  Dept  of  the  Navy 
Attn  OP  03EG 
Washington  DC  20350 


38 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  coiiection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  compieting  and  reviewing  the  coiiection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
coiiection  of  information,  inciuding  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1 21 5  Jefferson 
Davis  Highway,  Suite  1204,  Ariington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

December  1 997 


3.  REPORT  TYPE  AND  DATES  COVERED 


Progress,  from  1  October  1993-July  1997 


4.  TITLE  AND  SUBTITLE 


Light  Scattering  From  an  Optically  Active  Sphere  Into  a  Circular  Sensor 
Aperture 


6.  AUTHOR(S) 


J.  David  Pendleton,  David  L.  Rosen,  and  James  B.  Gillespie 


5.  FUNDING  NUMBERS 

DA  PR:  B53A 
PE:  P61102 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Laboratory 
Attn:  AMSRL-IS-EE 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1197 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Laboratory 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1197 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


ARL-TR-1474 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 

AMS  code:  611102.53A11 
ARL  PR:  7FEj60 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

Two  methods  are  developed  to  average  V  and  /  Stokes  parameters  over  a  circular  sensor  aperture 
collecting  light  scattered  from  an  optically  active  sphere.  One  method  uses  a  two-dimensional 
numerical  integration  of  Bohren's  theory  in  scalar  form,  which  explicitly  shows  the  connection  with 
Mueller  matrix  elements.  The  second  method  uses  expansions  of  vector  spherical  harmonics  in 
Bohren's  theory  to  integrate  over  apertures  of  any  size  without  convergence  checks  since  numerical 
integration  is  avoided.  Equations  simplifying  the  analytical  results  for  47t  scattering  are  obtained. 
Sample  computations  of  average  Stokes  parameters  are  performed  with  both  methods,  and  agreement 
to  six  decimal  places  is  obtained. 


14.  SUBJECT  TERMS 


Optical  activity,  circular  dichroism,  optical  rotation,  chiral,  polarization, 
Stokes  parameters,  Mueller  matrices,  Mie  theory,  biodetection,  biosensing 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 


Unclassified 


Unclassified 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

Unclassified 


15.  NUMBER  OF  PAGES 

44 


20.  LIMITATION  OF  ABSTRACT 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2-89) 
Prescribed  by  ANSI  Std.  Z39-18 
298-102 


39 


